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Abstract 

We investigate the physical properties of a two-component virialized protogalaxy com- 
prising a hot gas and cold clouds which are in pressure equilibrium under the assumptions that 
the protogalaxy is both spherical and homogeneous. Two conditions which we adopt are that 
the protogalaxy is in virial equilibrium and that the cooling time is less than the dynamical 
time. It is found that there exists an upper mass bound for such a two-component protogalaxy 
when the cloud mass is comparable to or greater than that of hot gas. We also consider the 
stability of the cold cloud against the Jeans instability and Kelvin-Helmholtz instability. 
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1. Introduction 

Recently, a remarkable series of observations has revealed many high-redshift galaxies, 
radio sources, and quasars to the redshifts greater than four, thus suggesting important impli- 
cations about galaxy formation and evolution (McCarthy 1993). Theoretically, from the 1980s 
many large simulations concerning galaxy formation in the context of large-scale structure 
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formation have been performed (e.g. Efstatiou, Silk 1983; Cen, Ostriker 1992 and references 
there in). Since there still remain many uncertain processes in galaxy formation, it could be 
valuable to clarify the fundamental nature of galaxy formation by a simple analytic treatment. 

A few decades ago several attempts were proposed (Ostriker 1974, talk at 7th Texas Con- 
ference; Gott, Thuan 1976; Rees, Ostriker 1977; Silk 1977; Rees 1978). In these studies, it was 
proposed that a protogalaxy consists of one component homogeneous gas in a virial equilib- 
rium state; it was assumed that radiative cooling controls the fate of the protogalaxy, with the 
ratio of cooling time to the collapse time being a controlling parameter. Their conclusions are 
that if free-free cooling dominates over bound-free cooling, the upper bound to the radius of a 
protogalaxy can be uniquely expressed by only fundamental constants, and that if bound-free 
cooling predominates, its mass is uniquely expressed as well. 

On the othe hand, another model was proposed in which the protogalaxy comprises a hot 
gas and cold clouds (Fall, Rees 1985) because the hot gas component is thermally unstable to 
make cold clouds with mass around globular clusters. Regarding such a two-component model 
of a protogalaxy, Ikeuchi and Norman (1991) examined the equilibrium structure, and obtained 
the expressions of its physical quantities by using fundamental constants and several additional 
parameters. Recently, Mathews and Schramn (1993) and Lee et al. (1994) have studied the 
star-formation history of a two-component galactic halo within the expanding background 
universe. They emphasized that in order to explain the age discrepancy in the galaxy two 
different types of star-formation modes are required. One is such an efficient star-formation as 
induced by cloud-cloud collisions in the halo; the other is the quiescent star-formation in the 
disk. 

In this paper we consider the physical properties of a two-component virialized protogalaxy 
by a simple analysis, and give an upper mass bound for it when the cloud mass dominates over 
the hot gas mass. We also consider the stabihty of cold clouds against the Jeans instability 
and Kelvin-Helmholtz instability. 
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2. An Upper Mass Bound 



2.1. Basic Equations and Assumptions 

Following the discussion by Ikeuchi and Norman (1991), we assume for simplicity that 
this two-component protogalaxy is spherical and homogeneous. Suppose that Nc cold clouds 
of each radius Rc, mass M^, and temperature are embedded in a hot gas of radius i?h, 
mass Mh, and temperature T^. To determine the physical properties of this two-component 
protogalaxy, seven equations, or reasonable assumptions for seven quantities in the above, are 
needed. They are as follows: 

(i) Virial equilibrium, 

al = 0.6 ^ „ ^ = 3 -, (1) 

Rh fJ'hrrip 

where (Jv(Th),/Xh, and rup are the virial velocity of protogalaxy, the mean molecular weight 
of hot ambient gas, and the proton mass, respectively. For a primordial wholly ionized gas, 
/ih = 0.609. We do not explicitly include the gravity of dark matter in this paper. 

(ii) Equality of the free-fall time and cooling time (rfj = Tcooi), 
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32 ^ G(Mh + iVeMc) • MhCAA(rh) 

where = An/S is a volume element, A(Th) is the cooling function of a primordial gas 
in units of erg s~^ cm'^ (Binney, Tremaine 1987) and c\ = 0.821 is a correction factor to 
the fitting formulae for the cooling function. Strictly speaking, the protogalaxy will collapse 
when the condition tq > Tcqoi is satisfied. We should thus note that equation gives the 
maximum mass/minimum radius of a collapsing protogalaxy. If the luminous mass of the 
galaxy is condensed in the dark halo (White, Rees 1978), we must consider a vast amount 
of dark matter. A zeroth-order incorporation of this effect can be achieved by changing the 
gravitational constant G to G[H-Md/Mb] in equations (|I|) and (0), where and Mb are the 
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masses of the dark matter and baryonic matter, respectively. 

(iii) Pressure equilibrium between clouds and hot ambient gas, 

= ^Te, (3) 

where /ic is the mean molecular weight of the gas in cold clouds of primordial abundance, 
/^c =1.23. 

(iv) Energy balance condition, 

nlclA{T^)c,{Rl-N,Rl) = Jk^^Rla.rjiN.M.c'). (4) 

Cv-Kh 

Here, we assume that the energy loss from the hot region is supplied by the energy liberation 
by the cloud-cloud collisions, and rj denotes the fraction of the liberated energy to the cloud 
rest mass. In the following, we suppose that rj ^ 10~^ — 10~^. These values are the observed 
ones in molecular clouds, considering that the mass fraction of energy liberation from newly 
born stars is ~ 10"'^, and the star-formation efficiency from clouds is in the order of 10~^ — 
10~^. We suppose that this cloud-cloud collision induces efficient star-formation during an 
early stage of protogalactic evolution. For our model to be consistent, the cold clouds must be 
stable (discussed in section 3). 

Conditions (i) — (iv) are the physical assumptions for a two-phase system. To determine 
seven quantities, we need to specify additional three conditions, though we do not have appro- 
priate ones at present. We regard the three quantities T^, T^. and as being free parameters, 
where is the volume-filling factor of clouds, 

F. = ^. (5) 

Using the characteristics of the cooling function of primordial gas, we suppose that the hot 
gas temperature Th is > in which the hot gas is thermally stable, and that the cold 

clouds temperature Tc is lO^K, in which the primordial cooling function has a sharp cut-off. 
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because of recombination to hydrogen atoms. A more detailed discussion concernig the cloud 
temperature is given by Kang et al. (1990). 

By using the above four equations (|l|) — (§) with three free parameters -Fv,Tc, and Th, we 
can obtain the equilibrium masses and radii of the hot region and cold clouds as follows: 

- "-^^ ^ cm^" 

n - 0440^^ x^MM!^ (H) 

^6 m0.5A(rh)ri-5(r/c2)3 



where 



/ = (l-F.)+2.02F.p^. (10) 
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2.2. Results 

In figure 1, we show the relation between the total mass Mtotai(= Mh + NcMc) and the 
radius -Rh of a two-phase protogalaxy for = 10^ K, rj = 10~®, and five cases of F^. For a 
comparison, we plot the case of a one-component protogalaxy by a dashed line. As is well- 
known, in the case of a one-component protogalaxy (no clouds) the upper bound of the radius 
is obtained at the free-free cooling region. On the other hand, in the case of a two-component 
protogalaxy the total mass is limited lower with increasing the cloud fraction. For example, 
the total mass is bounded to 10^'^ Mq for F^ = 10~^. For this result, we can give a clear physical 
interpretation. The cooling time is determined by the hot gas, while the dynamical time is 
controlled by clouds for a cloud-dominant galaxy. That is, in a cloud-dominant galaxy the 
protogalaxy, itself, cannot collapse, because of an inefficient energy loss from the hot region 
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[see equation (^]. The upper mass bound which we derived here is about two orders larger 
than the typical mass of a galaxy, ~ IO^^Mq. We discuss this point in section 4. 



3. Stability of Clouds 

In this two-component model, we assume that cold clouds collide with each other and 
liberate energy to the hot ambient medium. In order for this assumption to be satisfied, the 
cold clouds embedded in the hot medium must be dynamically stable. Therefore, we pose two 
stability conditions concerning the clouds. 

3.1. Gravitational Instability 

First, clouds must be gravitationally stable, 

M,<Me,jeans,p, (H) 

where the Jeans mass of a pressure-confined cloud, Mc^jeans.p, is given by (Spitzer 1978), 
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Here, Cc and are the sound speed in the cloud and the pressure of the ambient hot medium, 
respectively. We apply (p!2| ) to the cold cloud. 



c,Jeans,p y ^^O^K^^ Th ' 4oOkpc^ V Mh ° ^ ' 

This mass is roughly comparable to the mass of a globular cluster (Fall, Rees 1985). Of course, 
if we take Tc = 10^ K, due to II2 cooling ( see Kang et al. 1990), this situation greatly changes. 
In this case, according to equation (|13D, the Jeans mass of a cold cloud is smaller by four orders 
of magnitude than that for the case Tc = 10^ K. 

3.2. Kelvin- Helmholtz Instability 

Second, clouds must be stable against the Kelvin-Helmholtz instability (Chandrasekhar 
1961; Murray et al. 1992). This instability occurs when two media having different densities 



are in relative motion. In the present situation, clouds with a density Pc^Ph niove in the hot 
medium with random velocity a^. This instability occur when the momentum from a cloud to 
hot gas is transferred (cf. Chandrasekhar 1961), i.e., 

2pcRc > PhRh. (14) 

In another description, it reduces to 

1 p T 

Rc > -Rc,crit = - — Rh = 0.248 — i?h 

2 Pc Th 



Tc ,,10" K, / R^ 



= 2.48xlO^(^^)(^^-^)J ° pc. (15) 
404K^^ Th ^yiOOkpc ^ ^ ' 

Since in this paper we assume that the cold clouds are confined by the pressure of the hot 

ambient medium, the stabilizing mechanism against the Kelvin-Helmholz instability, such as 

the self-gravity and magnetic field of the cold cloud, can be neglected, and the above simple 

treatment for the clouds is not sufficient. However, this is the most generous criterion. We 

call equations ( p!TD and (|15]) the stability conditions, and Mc,jeans,p and -Rc.crit the critical mass 

and radius for the cold clouds, respectively. 

3. 3. Results 

In figure 2a, we show the mass of a cold cloud as a function of the hot gas temperature 
Th for Tc = 10^ K, = 10~^ and five cases of F^. The solid and dashed lines show M^. and 
Mc^jeans,p5 which are calculated from equations (|) and (|13]), respectively. In figure 2b, we show 



the radius of the cold cloud as a function of the hot gas temperature for the same parameters 
as in figure 2a. The solid and dashed lines show Rc and -Rc,crit, which are calculated from 
equations (|[) and (|1^), respectively. For different Fv, Mc,jeans,p and -Rc,crit change only slightly 
while Mc sensitively changes. The gravitational stability condition gives the upper limit of 
the cloud mass (or Th), and the Kelvin-Helmholtz stability condition gives the lower limit of 
the cloud radius (or Th). As can be seen in figure 2, there are only small allowed regions of 
parameters for each F^. Table 1 summerizes the allowed ranges of Th for each model, and 



table 2 gives the physical quantities of a two-component protogalaxy when the median value 
of Th from table 1 is taken for each model. In our model, the cold clouds of mass ~10^Mq are 
confined by the thermal pressure of a hot ambient gas of mass ^lO^^M©, irrespective of F^. 

4. Summary and Discussion 

We considered the physical properties of a two-component protogalaxy by using a simple 
analytic treatment, and determined them considering the stability of the cold clouds. In the 
cloud-dominant protogalaxy, the mass of a protogalaxy is almost constant in the free-free 
cooling region above ~10^ K. This means that a protogalaxy heavier than ^lO^^M© cannot 
collapse. In contrast, the mass of an one-component protogalaxy is not limited. The stability 
condition for cold clouds against the Jeans instability and Kelvin -Helmholtz instability poses 
stringent limitations on the physical quantities of the two-component protogalaxy. 

In our simple model, several important factors are neglected concerning the formation 
of galaxies, such as the abundance and distribution of dark matter (White, Rees 1978), the 
cosmological initial conditions, the effect of shock heating (Thoul, Weinberg 1995), the angular 
momentum and the deviation from spherical symmetry of a protogalaxy (De Araujo, Ophcr 
1994). To explain the typical mass of a galaxy (~ IO^^Mq), quantitatively, we must take into 
account the above-mentioned factors. The upper mass bound which we present here may give 
one possibility for the origin of the typical galaxy mass. 

According to the present viewpoint of a two-component protogalaxy, we speculate on the 
star-formation history as follows. In the protogalaxy with radius i?h ~ 100 kpc and a total 
mass of Mtotai ~ lO^^M©, the star-formation in an extended halo is triggered by cloud-cloud 
collisions and/or the shock compression of clouds by the jet suggested in high-redshift radio 
galaxies (McCarthy 1993; Norman, Ikeuchi 1994, private communication). These Population 
II stars distribute over the extended halo, and their low mass ends may be candidates for 
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MACHOs. Meanwhile, the hot gas is cooled and contracts to a disk or spheroid, depending 
upon its angular momentum. Prom this gas component. Population I stars are formed. 

In the succeeding work we consider the evolution of this two-component protogalaxy. Es- 
pecially, once the protogalaxy begins to collapse, the hot gas phase is cooled and contracts, 
while the cold clouds can not contract, because their random energy can not be extracted. 
That is, the hot gas and cold clouds segregate. At this stage the stability of the cold clouds 
changes, because of the lack of hot ambient gas. Including these processes, we examine the 
evolution of two-component protogalaxies. 

The authors would like to thank Professors M.Sasaki and N.Gouda for valuable discus- 
sions, and K.M. is grateful to Dr.M.Shibata for helpful suggestions and a critical reading of 
the manuscript. 
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Table 1. Allowed values for four parameters with respect 
to the Jeans and Kelvin- Helmholtz instability. 





Model logTc 


log 7^(10- 


-«) Fv(10-4) 


log Th(K) (allowed) 




1 




3.0 


6.38— 


6.41 




2 


10 


2.0 


6.73— 


-6.78 




3 

4 




1.0 


7.33— 


-7.40 




4 




10.0 


6.34— 


-6.42 




5 


1 


8.0 


6.53— 


6.64 




6 




5.0 


6.94— 


-7.12 




Table 2. Physical quantities of a two- phase protogalaxy 




for the median value of Th in 


table 1. 




Model logTh(K) 


log i?h(pc) 


logMh(MQ) 


logi?c(pc) 


log Me (Mo) 


1 


6.40 


5.39 


12.93 


2.41 


6.68 


2 


6.74 


5.29 


13.14 


1.95 


6.18 


3 


7.36 


5.13 


13.51 


1.18 


5.37 


4 


6.40 


5.28 


12.70 


2.34 


6.59 


5 


6.60 


5.19 


12.77 


2.06 


6.28 


6 


7.00 


5.03 


12.92 


1.49 


5.60 
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Figure Captions 

Fig. 1 

Relation between the total mass and radius of a two-component protogalaxy for Tc = 10'' K, 
rj — 10~^ and five cases of Fy. For a comparison, the case without clouds is also plotted by a 
dashed line. 
Fig. 2 

Cloud mass (a) and radius (b) as a function of Th for Tc = lO'' K, rj = 10~^ and five cases of 
Fy. The Jeans upper mass limit and the Kelvin-Helmholtz lower radius limit are illustrated 
by dashed lines for the same parameters. 



12 




12 12.5 13 13.5 14 14.5 

logMtotQiCMo) 




4 4.5 5 5.5 6 6.5 7 7.5 8 



logT, (K) 



